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Abstract 

An interfacial approximation of the streamer stage in the evolution 
of sparks and hghtning can be written as a Laplacian growth model 
regularized by a 'kinetic undercooling' boundary condition. We study 
the linear stability of uniformly translating circles that solve the prob- 
lem in two dimensions. In a space of smooth perturbations of the 
circular shape, the stability operator is found to have a pure point 
spectrum. Except for the eigenvalue Aq = for infinitesimal transla- 
tions, all eigenvalues are shown to have negative real part. Therefore 
perturbations decay exponentially in time. Wc calculate the spectrum 
through a combination of asymptotic and series evaluation. In the 
limit of vanishing rcgularization parameter, all eigenvalues are found 
to approach zero in a singular fashion, and this asymptotic behavior is 
worked out in detail. A consideration of the eigenfunctions indicates 
that a strong intermediate growth may occur for generic initial pertur- 
bations. Both the linear and the nonlinear initial value problem are 
considered in a second paper. 

PACS: 47.54.-r 

Keywords: moving boundary, kinetic undercooling rcgularization, lin- 
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1 Introduction 



The motion of interfaces in a Laplacian field is of general interest and has 
been a subject of intense study over many years (see for instance the reviews 
[U [2]). Such problems arise in many physical contexts, such as viscous 
fingering in multi-phase fluid flowPllllIlll[7l[S[Sll0l[IIl[l2l[l3lllllll5l 
I16j . dendritic crystal growth in the quasi-steady small Peclet number limit 
[13 H 13 [22] , void electromigration [Ml [23 [23 EZ] and a host of other 
phenomena such as the growth of biological systems like bacterial colonies 
or corals [2H]- 

More recently, a similar mathematical problem has been discovered in 
the study of 'streamers' [29l[30l[3ll[321[33l[3ll[35l[36] which occur during the 
initial stage of electric breakdown and play an important role both in the 
natural phenomena of sparks and lightning as well as in numerous technical 
applications [35]. Streamers are weakly ionized bodies growing into some 
nonionized medium due to an externally applied electric field. This field is 
so strong that the drifting electrons very efficiently create additional electron 
ion pairs by impact ionization, and the nonlinear coupling between ionized 
body and field further increases this effect. 

Models for negative streamers in simple gases like nitrogen or argon are 
based on a set of partial differential equations for the densities of electrons 
and of positive ions coupled to the electric field [301 [311 [321 [331 [311 [35] . 
Analysis and numerical solutions of these equations reveal that in the front 
part of the streamer, a thin surface charge layer develops where the electron 
density strongly exceeds the ion density. Therefore the electric field E varies 
strongly when crossing this layer. Right before the layer, it is enhanced, but 
in the interior of the streamer, it is screened to such a low level that impact 
ionization is suppressed and the electron current transporting charge from 
the interior to the surface charge layer is small. Consequently we may take 
the interior as being essentially passive, and the growth of the streamer is 
governed by the surface charge layer which is driven by the strong local field. 

If the external field is very strong, the thickness I of the surface charge 
layer can become small compared to the typical diameter 2i? of the streamer 
|36| . This suggests modeling this layer as an interface separating the ionized 
interior from the nonionized exterior region. In this model the variation of 
the potential ip across the surface charge layer is replaced by a discontinuity 
on the interface. Since the interior is considered passive, only the limiting 
value (/J"*" reached by approaching the interface from the outside is relevant for 
the dynamical evolution, and analysis of results of the PDE- model suggests 
[361 [371 [38] that with an appropriate gauge, (^"'" is coupled to the limiting 
value E"*" of the electric field by the boundary condition 

(^+ = -^n • E+ . (1) 

Here n is the outward normal on the interface, I is the regularization length 
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corresponding to the interface thickness, and E'^ = — V^j^, where the 
again indicates the Hmit of approaching the interface from the outside. In 
the context of dendritic crystal growth, this boundary condition is equivalent 
to including the kinetic undercooling effec10 while excluding the usual Gibbs- 
Thompson surface energy correction to the melting temperature. 

As the motion of the interface is caused by the drift of the electrons in 
the local electric field v = — E, the interface moves with normal velocity 

Vn = -n- E+, (2) 

and outside the streamer the potential obeys the Laplace equation 

Aip = 0. (3) 

We discuss the problem defined by Eqs. ([I])-® in infinite two-dimensional 
space, with the electric field becoming constant, E = — V(/> — > Eqo, far 
from the streamer; such a condition is realized frequently in atmospheric 
discharges, e.g., inside thunderclouds. This far-field condition on ^(p is 
different from the usual source/sink condition in viscous fingering in the 
absence of side- walls, or undercooling specification in the quasi-steady low 
Peclet number crystal growth problem which have been extensively studied. 
However, moving bubbles and fingers in a long Hele-Shaw channel are in- 
deed subjected to this type of far-field condition. (For the discussion of the 
streamer problem equivalent to the Saffman- Taylor finger, we refer to |4Uj.) 
Further, in the crystal growth or directional solidification problem, a sys- 
tematic inner-outer analysis for small Peclet number |17[ [21] shows that this 
is an appropriate condition at co for the "inner" problem. 

A simple steady solution to the streamer equations given above is a circle 
translating with constant velocity determined by Eoo [371 138] . Though such 
circles differ from proper streamers, which are growing channels of ionized 
matter [351 HOj WL\ I42j . their front half closely resembles the head of the 
streamer where the growth takes place. It is therefore a question of physical 
interest whether or not translating circular solutions are stable to small 
perturbations and this is the subject of the present investigation. 

The relevance of this analysis for more realistic streamer shapes is sup- 
ported by results found in another physical context. Steadily translating 
circles also arise in viscous fingering in a Hele-Shaw cell when surface ten- 
sion is included (instead of kinetic undercooling) in the limit when the bubble 
is small compared to the cell-dimensions [3]. The linear stability of these 
bubbles, including larger non-circular steadily translating bubbles, has been 
studied before both for one and two fiuids \10\ [T3] and the results largely 

^Under most natural conditions of crystal growth, kinetic undercooling is important in 
a limit when the Peclet number is not small enough to justify a Laplacian field approxima- 
tion; nonetheless, there have been some studies of steady Laplacian crystal growth with 
kinetic undercooling effects only [39] . 
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mimic those obtained for a finger, though the latter calculations are math- 
ematically much more involved. 

It has been known for a while that in the absence of any regularization, 
such as surface tension or kinetic undercooling, the initial value problem 
in a Laplacian field is ill-posed [9] in any norm that is physically relevant 
to describing interfacial features. This is reflected in the instability of any 
steady shape, when the growth rates increase with the wave numbers of the 
disturbances. Ill-posedness makes idealized model predictions sometimes 
physically irrelevant (see [TB] for a thorough discussion) and regularization 
becomes essential. 

Considering a planar front, one finds that regularization does not remove 
the instability against fluctuations of small wave number. However, for large 
wave numbers, linear stability analysis exhibits a basic difference between 
surface tension and kinetic undercooling. All large wave number components 
of a disturbance decay with surface tension regularization, while for kinetic 
undercooling the growth rate saturates to a constant that scales as [16j : 
for streamers, such a saturating dispersion relation is derived and discussed 
in [13111]. 

For curved fronts, one can pose the question: how does curvature sta- 
bilize, if at all, a disturbance whose wavenumber is in the unstable regime 
for a flat interface, either with surface tension or with kinetic undercool- 
ing regularization? With surface tension regularization, some answers are 
available in the existing literature. Arguments have been presented [H [6] 
that suggest that a localized wave packet with wave numbers in the unstable 
regimqj advects along the front as it grows; once it reaches the side of the 
front where the local normal velocity is zero, the disturbance stops growing. 
If the steady shape is closed, as it is for a circle, the continued advection 
of the wave-packet towards the receding parts of the interface will cause 
the disturbance to decay eventualljU. If regularization is small, there is a 
large transient growth. Unless the disturbance amplitude is smaller than a 
threshold that shrinks to zero with regularization, the transient exponential 
growth causes the interface to enter a nonlinear regime that can destabilize 
the steady front, even when it is predicted to be linearly stable. Analysis 
of approximate equations, supported by numerical calculations of the full 
equations support the above scenario. Similar stabilization should occur 
for the kinetic undercooling boundary condition as well, though we are not 
aware of any explicit study affirming this expectation. 

^Localized disturbances refer to those with wavelengths far smaller than the typical 
radius of curvature of the steady shape. These can be unstable only if the regularization 
parameter is sufficiently small. 

^Surface tension causes localized disturbances to decay as they advect to the sides 
even when an interface is not closed but becomes parallel to the direction of motion as 
is the case for a finger in a Hele-Shaw cell. However, no decay is expected for kinetic 
regularization. This is where a closed interface is different. 
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Note, however, that stabihzation of locahzed wave packets does not rule 
out instabihty to long ranged disturbances. A formal asympotic study for 
small nonzero surface tension [121 [Tl] as well as numerical studies [7j reveal 
that surface tension stabilizes precisely one branch of steady solutions for 
fingers and bubbles in a Hele-Shaw cell. Similar results follow for a needle 
crystal [18] though in the latter case, convective instability of wave packets 
caused by significant normal speed along the parabolic front is believed to 
cause dendritic structures [I]. These conclusions have been challenged at 
times by alternate scenarios (see for instance |19j ) that are based on for- 
mal calculations, but with different implicit assumptions. Such controversy 
affirms the need for more rigorous mathematical studies of the stability prob- 
lem, even if it is for relatively simple shapes such as the circle in the present 
study. 

For the kinetic undercooling boundary condition, relying merely on a 
numerical study to understand the long time behavior is fraught with dif- 
ficulties. One finds a collapsing spatial scale for large time at the rear of 
the circle. Analytically, this is found for e = ^/i? = lin [37^ 138 j^ . As 
will be argued in the present and the companion paper, the occurrence of 
this collapsing scale is a general feature for any e > 0. This means that 
as i — > oo, one must resolve progressively finer scales near the back of the 
bubble. Further, calculations for small e require resolving a large number of 
transiently growing modes. All this underscores the need of some progress 
on the analytical side. 

The present paper, which is part I of a two-paper sequence, is devoted 
to the spectral properties of the linear stability operator, associated with 
infinitesimal perturbations of a circle. In part II [35], we will consider the 
initial value problem, presenting analytical and numerical results on the 
evolution of both infinitesimal and finite perturbations. 

The present paper is organized as follows. In Section we reformulate 
the problem defined by Eqs. ([I])-© by standard conformal mapping, and we 
present the PDE governing the time evolution of infinitesimal perturbations 
of the circle. This material has been presented before in [37^ [38], where 
also the general solution of the PDE in the case e = 1 has been discussed 
in detail. The explicit solution found for e = 1 shows that outside any 
fixed neighborhood of the rear of the bubble, the long-term behavior of 
infinitesimal perturbations is described by X^^o ^"^"^/^An ■, where A„ is the 
nth eigenvalue (ordered according to absolute value) of the linear stability 
operator and /3a„ is the corresponding eigenfunction. 

We then study this eigenvalue value problem for arbitrary e > 0. We 
show in SectionOthat the linear stability operator, defined in an appropriate 
space of analytic functions, has a pure point spectrum. In Section S] it is 

*We recall that i? is a measure of the size of the streamer. The precise definition is 
given in Eq. ^ below. 



5 



proven that there are no discrete eigenvalues with non- negative real part, 
except A = that corresponds to the trivial translation mode. A set of 
discrete, purely negative eigenvalues is calculated in Section [5] as a function 
of e; they smoothly extend the results found previously for e = 1. The results 
suggests that as e — > 0, the spectrum degenerates to the trivial translation 
mode and this limit is discussed in detail in Section [H Section [7] contains a 
discussion of the eigenfunctions belonging to these eigenvalues, and Section [H] 
contains the conclusions. Some part of our analysis exploits general results 
on the asymptotic behavior of the coefficients of Taylor expansions. These 
results are presented in an appendix. 

2 Reformulation by conformal mapping 

In this section, we collect results and notations from (STJ [38] that will be 
used in later sections. 

2.1 Problem formulation and rescaling 

We consider a compact ionized domain D in the (x, y)-plane. We assume 
that the net charge on the domain vanishes (i.e., it contains the same number 
of electrons and positive ions). The domain moves in an external field that 
far from the domain asymptotically approaches 

Eoo = -|Eoo|x. (4) 

Here x is the unit vector in x-direction, and |Eoo| sets the scale of E and 
thus of the potential (p. As length scale we take 




where \T>\ is the area of T>, which is known to be conserved. This follows 
from the charge neutrality of the streamer since 

0= dx dyV ■'E= ds n{s) • E = - / ds Vn, (6) 
Jv JdV JdV 

where in the last step we inserted Eq. ^ for the normal velocity of the 
boundary. Since 

/ dsvn = dt\V\, (7) 

JdV 

the area is conserved, irrespective of the precise charge distribution in the 
interioiH. Also introducing the time scale -R/|£'oo|, we rescale the basic 

^We remark that the argument is straight forward to generalize to three spatial dimen- 
sions. Therefore the volume of a charge neutral object with surface velocity v oc E"'" in 
three spatial dimensions is conserved as well. 
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equations to the dimensionless form 



Ayp = 0, {x,y)(^V 
Vn = rv - (Vv7)+ 
= en-(V(^)+. 



(8) 
(9) 
(10) 



The only remaining parameter in the rescaled problem is 

e = l/R. (11) 
The boundary condition at infinity after rescaling takes the form 

(/9 — > X + const for yaP^+ip — > oo . (12) 

2.2 Conformal mapping 

We now identify the physical (x, y)-plane with the closed complex plane 
z = X + iy^ and we introduce a conformal map f(ij,t) that maps the unit 
disk in the w-plane to the complement of D in the z-plane, with a; = 
being mapped on z = oo, 



z = f{u,t) = ^^ + f{u;,t), a„i(t)>0. 



UJ 



We further define a complex potential <^(w,i) obeying 
Re[$(w, t)] = (p{f{uj, t)) for uj gU^ . 
The boundary condition (fT2]) and the Laplace equation ([8]) 



UJ 



(13) 
(14) 

enforce the form 
(15) 



with being holomorphic for uj G U^j. 

The two boundary conditions ([9]), (fTOl) take the form 



\d^f \ Re[^>] 



Re 



for UJ G dUrc 



d.fW 

eRe[a;9a;^>] for w G S^o; , 



(16) 
(17) 



which completes the reformulation of the moving boundary problem (l8l)- ()12p 
by conformal mapping. 

We will restrict the analysis here to initial conditions /(w, 0) holomorphic 
in some domain IAq D lA^. In part II [l5] of this paper sequence, we will give 
evidence that analyticity on lA^^j is preserved in time, though the distance 
of the domain of analyticity hl[ to dU^ shrinks with time. The streamer 
boundary 9D, which is the image of boundary dU^ under f(io,t), will turn 
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out to be analytic and therefore smooth. Similar analytic representations 
exist for the entire class of 2-D Laplacian growth, with details depending 
on the type of boundary condition, geometry and asymptotic conditions at 
infinity. For the classic viscous fingering problem, Polubarinova-Kochina 
[21j and Galin [20] use a representation that coincides with the one given 
above in the unregularized case e = 0. 

2.3 Linear perturbation of moving circles 

It is easily seen that equations ()16p . (I17p allow for the simple solution 



$(0)cU,t = ---—CO, 

LO 1 + e 

which in physical space describes circles of radius 1 moving with constant 
velocity 2/ (1+e) in x direction. (We recall that the radius was scaled to unity 
in Section [2m ) We note that relaxing the analyticity conditions on f{u},t) 
on = 1, one can obtain another set of uniformly translating solutions, as 
recently discovered jH]. The present paper is restricted to perturbations of 
the steady circle that retain the imposed analyticity of the streamer shapes, 
and hence analyticity of / (as well as $) on |u;| = 1. 

As the area is conserved (as shown in subsection l2.ip . the residue a_i = 1 
does not change to linear order in the perturbation. We therefore can use 
the ansatz 

/(u;,t) = f(^\u;,t) + rj (3iu;,t), 

(19) 

where 77 is a small parameter, and P{Lo,t), x(<-^,i) are holomorphic vuli^. A 
first order expansion of Eqs. ()16|) . ()17|) in 77 yields the following boundary 
conditions for the analytic functions f3{uj,t) and x(a;,t) on = 1: 

Re[ujdrf3 — oodi^fi] = Re[—Lodi^x\ , 

(20) 

^Re[{u; + ^)L0^d^f3] = Re[eLod^x + x] , 



r = ^^t. (21) 



where we rescaled time as 



Since the left and right sides of each of the two equations in (j20p are real 
parts of analytic functions and each is assumed a priori continuous upto the 
boundary, they can differ everywhere in lo by at most an imaginary constant. 
Evaluation at w = shows this constant to be zero for the first of the two 
equations. Elimination of x results in the linear PDE: 



CeP = (22) 
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with the operator 

^e = ^ (w^ - I) ijd^ + ed^cudr + dr - . (23) 

We note that is of similar structure as the operator resulting from a 
linear stability analysis of translating circles in the context of void electro- 
migration [24:\ [26] . The main difference here is the occurrence of the mixed 
derivative dojUjdr- 



2.4 Formulation of the eigenvalue problem 

To motivate our formulation of the eigenvalue problem, we note some re- 
sults on the temporal evolution of infinitesimal perturbations. In [38], the 
equation C^P = was solved as an initial value problem for the special 
value e = 1. It was found that any initial perturbation /3(a;,0) holomorphic 
in U' D Z^i^ for T — > oo is exponentially convergent to some constant. This 
results from the expansion 



P{co,T)=J29n(3^^!{io)e'"\ (24) 



n=0 

with 



Xn = -n, n e No, for e = 1. (25) 
The coefficients gn and the eigenfunction^ 



/t'H=/ ^ (^1 (26) 



X dx f X — 1 
Uj'^ \x + \ 

are determined by an expansion of {2+ujd^S)(i{'jJ-, 0) in powers of (1 — w)/(l + 
cj). For n > the eigenfunctions (f26l) are singular at w = — 1, though [3{uj, r) 
is not. The expansion ()24p is convergent in a domain 7)^ expanding in time 
that eventually includes every point in U^^ \ {—1}. For large r the region 
where the expansion is invalid, shrinks to lij = — 1 exponentially. This region 
is measured by the new scale 771(0;, r) = (1 + uj)e^ , and the expansion ([^^ 
is valid if r]i is large. For r/i < 0(1) the perturbation for r — > 00 behaves 
as I3{lo,t) Fq{i]i) + 0{e~'^) where Fq is some analytic function of its 
argument, depending on P{uj,0). 

For an analytic initial condition on U' , with a lone branch point singu- 
larity ujs in \uj\ > 1 not on the positive real axis, the emergence of this new 
scale near a; = — 1 can be related to the approach of this complex singularity 
towards —1 exponentially in r for large r. Asymptotic arguments that will 



®Note that the exponent — A„ in (|26p is correct while +A„ in Eq. (4.20) in [38] is a 
typo. 
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be presented in part II [l5] suggest that this behavior is generic for all e > 0. 
The analysis is based on the linear the evolution equations for 6^, where 



k=0 

For A; ^ e"^, we find the asymptotic relation 



6fc~(-lffc-"/i(T)exp [-kf{T)] 

where 

/(r) = log 



l + Ce 



l-Ce~ 



with C ^ ~^ 



UJ. - 1 



For LOs ^ (1,00), /(r) stays finite and approaches exponentially in r for 
large r. If lJs G (1,oo), /(r) increases monotonically to 00 for r G (0,rc) 
where e~'^'' = 1/C. For r > Tc, /{t) decreases monotonically and approaches 
exponentially in r as r — > 00. In either case, from the known relation 
between Taylor series coefficients and the location of the closest singularity of 
an analytic function (see the appendix), it follows that /(r) ~ e~'^ as r — > 00 
implies that /3 has a singularity approaching lj = —1 exponentially in r for 
large r. This feature is retained for any other isolated initial singularities as 
well, though /c~° is replaced by a more complicated dependence in k. Since 
the problem is linear, the evolution of a distribution of initial singularities 
can be understood from the linear superposition principle. 

This suggests that for any e > 0, as for e = 1, P{u!,t) has a collapsing 
scale (1 + uj)e'^ , and an expansion of the type (1240 cannot be valid in this 
neighborhood of u; = —1. 

Thus, in seeking an eigenfunction by substituting 

/3(^,r)=/3^)He"^ (27) 

into (f22]) . (f23]) . it is appropriate to allow Z?^*^^ to be singular at lo = — 1. 
Indeed, substituting the form (p7|) reduces Eqs. ([22]) . ([23]) to the eigenvalue 
problem 

L(e,A) /3j^(^) = 0, (28) 
L{e, A) = Ql + - 1) 5. + A(l + e + eujd^). (29) 

Evidently this ODE has three regular singular points, namely u; = and 
a; = ±1. The independent solutions at these points for e > are in leading 
order 

UJ 

[ITUJ] 



/?rH ~| J^2/e foru;^0, (30) 

-I nt^'iLx forc^il, (31) 
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We require the eigenfunctions (3^ to be solutions of ()28|) that are analytic in 
a; = and lo = 1. This is also the natural choice from a physical point of 
view since it is the right half of the circle, Re [a;] > 0, that corresponds to the 
physically interesting tip of the streamer. In general, eigenfunctions cannot 
be expected to be regular at all three points. Starting with a function regular 
at u; = 0, we cannot generally require regularity at both points uj = ±1 by 
adjusting the single parameter A. As shown in subsection 14. 3t the only 
eigenfunction regular at all three points is the trivial translation mode 



As noted above, an operator similar to P3p occurs in the problem of 
void electromigration, see section 4.1.3 in [26]. Again an eigenvalue analysis 
would yield a second order linear operator with three singular points at 
a; = and ±1 and therefore the eigenmodes in general cannot be regular 
at all three singular points. It is interesting to note that the authors 
conclude that their problem is unstable because the initial value problem for 
large time is singular at a; = — 1. In the current problem, the solution [l5] 
of the initial value problem is not singular at w = — 1; the singularity of the 
eigenfunctions does not reflect the true behavior of solution since, as has been 
pointed out earlier, there is an anomalous contracting scale e^{l+uj) near the 
back of the bubble. Whether or not there is an analogous contracting scale 
for the void electromigration problem [26] remains an interesting question. 
This anomalous scale shows up when the limiting processes lim^^_i and 
limT-_^-i„oo do not commute for the solution of the initial value problem. 

3 Discreteness of the spectrum 

We define A to be in the spectrum, if the linear operator L(e, A) does not 
have a bounded inverse in the class of functions / that are analytic in an 
arbitrary compact connected set V C ^'\{ — 1} that contains the whole 
line [0, 1] in its interior. A is in the discrete spectrum if ()28p has a nonzero 
solution p{^\uj) that is analytic in any such domain V. We now argue that 
if A is not in the discrete spectrum, then L(e, A) has a bounded inverse, i.e. 
there is only a discrete spectrum in this problem. 
To determine L~^, we solve the equation 



for a given h analytic in V, imposing the condition that also g is analytic in 
V. The solutions of the homogeneous equation L(e, A)/ = that are regular 
at u; = or u; = 1 will be denoted by fi{uj) or /2(w), respectively. It follows 
from Eqs. ()30p . (I3ip that these functions are determined uniquely up to a 
multiplicative constant. In the exceptional case where both independent 



Ao = 



/3q (cj) = const. 



(32) 




(33) 
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solutions are regular at = 1, A belongs to the discrete spectrum, see 
Section [5l A standard calculation shows that Eq. (j33p is solved by 

= r^^' G{uj,J) h{oj') + aM fi{Lo), (34) 

where 

'^^^^^'^= [/,(.)A(.')-/iM/2(-')], (35) 

and the coefficient ai[h] is a functional of h{uj'). C(A,e) does not vanish 
since otherwise the Wronskian fi d^f2 — f2 d^ifi vanishes identically and A 
is part of the discrete spectrum. It is easily seen that Eqs. (fMI) . (f35|l render 
g{uj) analytic in w = 0, and this condition eliminates any contribution of 
the form 02 [/i] f2{'-^)- 

Analyticity at u; = 1 is enforced by a proper choice of ai[/i]. To make 
the analysis explicit, in addition to /2(w), we introduce another solution to 
L[e, A]/ = by requiring 

h{u,) = {l-u:Yl^-^ h{^). (36) 

where fsiuj) is analytic at u; = 1. Using this form of fsiio), we exclude the 
case I — A G Z"^, that will be discussed later. Writing fi{u>) as 

/l(c^) =C2/2(c^)+C3/3(t^), (37) 

we find that G{uj,uj') from Eq. ((35]) takes the form 



G{UJ,LU ) 



(1+^')1A+A 



/2(.;)/3(a;') - ( h{u:)h{J) 



Evidently the first part in the square brackets for a; — > 1 yields a regular 
contribution to g{uj) from Eq. The contribution to J G h that is 

singular in = 1 has the form 

Jo 

where 

Hi-') = jYTW^' ^'^"'^ '^"'^ ^''^ 

is regular at c^;' = 1. If Re A — > —1, we can write 

-c^hiio) r dj {i-jf~^/^H{j) 

JO 

= -cM-) C dJ {l-Jf~^/'H{J) 
Jo 

-1 /I ,/xA-l/e 



+csh{u) J dJ [y^) ^(^')- (39) 
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The second part is regular at u; = 1 and the singular first part is canceled 
by the choice 

ai[h]= [ duj' {l-oj')^~^/'H{uj'). (40) 



We note that this result is valid also for A = ^ + n, n G N, where f^{oj) 
instead of being of the form (|36p shows a logarithmic singularity. 

If —n > Re A — i > — n — 1, n G N, we carry through n subtractions of 
H{u)') at io' = 1, defining 



n-l 



[H{J)\^ = H{J)-Y,H, {l-u^y, (41) 

j=0 

so that [H[u)')]^ ~ const (1 — w')". A short calculation shows that the 
singular part of J G h \s canceled by the choice 

aM = cL.' (1 - [H{u:')\^ + (^2) 

The expressions above clearly remain valid when ^ — Re A = n, except when 
i — A = n, a positive integer. 

When i - A = n is a positive integer, from well-known theory [U] for 
regular singular points, instead of (jSip . the solutions /i and /2 as defined 
earlier must have the following local representation near u = 1: 

/i(u;) = Ci(l-a;)'^Si(w)log(l-cj) + S2(c^), (43) 

/2(cu) = (l-u;)"i?i(w), (44) 

where Bi and B2 are analytic at a; = 1. If /i and /2 are independent, as 
they are when A is not in the discrete spectrum, then Ci(e, A) / 0. 
We now define 

H{u) = B2{u){l + w)"-2Aa;2A/i(a;), (45) 
while Hj is still defined by the expression (|41|) . It is also convenient to define 

Q{uj) = Bi{uj){l + uj)''-^/'u:'^''h{uj), (46) 

Note that each of H and Q are analytic at w = 1. Straight forward calcula- 
tion based on shows that the possibly singular part of g{uj) at u; = 1 is 
given by 

-^ln(l-u;)Si(.;)(l-a;)'^ 



+ fi{uj) i^i - j dJ 



l-UJ 
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The last term is analytic at w = 1. The singularity vanishes if we choose 

For any A for which C(A, e) ^ 0, using the explicit expression (I34p with 
ai determined from (j40p . ()42p or (j47p . whatever the case may be, we have 
in the domain V, 

Iblloo < C||/i||oo, (48) 

This conclusion follows from observing the properties of the integrand and 
noting that the Hj^ j = 1, ...n are bounded by some multiples of sup^^gy |/i(w)|, 
since they involve only a finite number of derivatives of /i at u; = 1. Since the 
boundary dV is at a finite distance from to = 1, the derivatives dijh\^_-^ by 
Cauchy's theorem are bounded by bj sup^gy I^C'^)! where bj is independent 
of h. Hence, we have showr0 that A is in the spectrum only if C(A, e) = 0, 
i.e. we can only have discrete spectrum in this problem. 

4 Absence of eigenvalues with positive real part 
and of purely imaginary eigenvalues 

4.1 Purely positive eigenvalues 

Real eigenvalues A > easily are excluded. Substituting into Eq. ([^5]) the 
power series 

OO 

/3j)(^) = (49) 

fc=0 

which converges for < 1 due to the location of the regular singular points, 
we find the recursion relation 

bk = 2A J^'^, fcfc-i + e 1^ forA:>2, (50) 

k{2 + ek) 2 + ek 

and 

bi = 2A i±^6o . (51) 
z -\- c 

We choose bo = 1 as initial value. 

For A > 0, evidently all bk are positive, and bk obeys the bound 

k — 2 
bk > e— — - bk-2 , 
2 + ek 



and therefore 



r {-) 

bk > —, V I., const. > 0. (52) 

r(i + i + |) 



^Note from definition, A is not in the spectrum if the resolvent L ^ is bounded. 
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For 1/e, this yields the lower bound 

bk > const k~^~^^^ , 

which shows that a sufficiently high derivative of I3^\uj) (j49p diverges for 
a; = 1, which contradicts the regularity requirement. 

4.2 Eigenvalues with positive real part 

To eliminate eigenvalues A = ^. + iv with > needs more refined argu- 
ments. We first derive an inequality replacing (I52p above. Motivated by 
(|52p . we rewrite the recursion relation (|50p in terms of 



Ck = ^^Vm^'^ ^ - ^^^^ 

This yields 

Ck = Xgk Cfc-i + Cfc-2, ^ > 3, (54) 

where 

^'^^ fc(2 + e^)r(|)r(i + i + |) "A^^UJJ- ^''^ 

We now multiply (j54p by c^_^ and take the real part. With the notation 

rfc = Re {cfc4_i} , (56) 

we get the relation 

rk = fJ- gk Ick-il"^ + rk^i, k>3. (57) 

Since ^ > and 

r2 - Keic2Ci| - 2^|A| (2+7)2 F(I) ' 

the Tfc form an increasing series of positive numbers bounded by 

rk>r2> 0. (58) 
The recursion relation (j57p formally is solved as 

k 

rk = tJ''^9j \cj-i\'^ + r2, k>3. (59) 
i=3 

Using now the relation 

|cfc|2 + |cfc-ip = 2rfc + \ck - Cfc_i|2 > 2rfc, 
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we find the bound 



^ > H2^gk\cj-i\ +r2. (60) 

We note that this bound is positive and increases monotonically. 

We now recall that by definition of the eigenfunctions, the only singu- 
larity in the complex plane is of the form 

(1+..)V^+V(^), 

where f{uj) is regular in \lo\ < 1, including uj = —1. A standard result on 
the relation of power series coefficients to the closest complex singularity 
(see the appendix) is, that the asymptotic behavior of of the expansion 
(HH} satisfies 

where -Boo is some constant. In view of (|53p . this implies 

\ck\ ~ const 2-1-^/' k-f" 

for /c — > oo which contradicts the bound (j60p . We thus conclude that there 
are no eigenvalues A with Re A > 0. 



4.3 Purely imaginary eigenvalues 

Now consider the possibility of a purely imaginary eigenvalue X = iv with 
z/ real. From the recursion relation, it is clear that 1^ = corresponds to 
the translation mode. So we only consider the case ^ 0. From complex 
conjugation symmetry, it is clear that if A = ii^ is an eigenvalue, so is 
A = —iu. Therefore, we may assume without any loss of generality that 
> 0. 
We introduce 

h = (61) 
The recursion relation (j5U|) takes the form 

h = 2. JL±^ _ , Izl bo = 1, = 0, (62) 

which shows that is real for any k. Thus the function 

oo 

P{u;) = ^bk6j'' = (3['\u), d; = iu;, (63) 
fe=i 

is real for real a), and the reflection principle guarantees 
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Thus /3(cj) either is singular both at u) = iti, corresponding to w = ^pli oi' 
is entire. When A = ii^ is some eigenvalue, we cannot have a singularity 
at uj = —i and so must be entire. We will now show that this is 

impossible. 

First, note from the recursion relation that (3{uj) cannot be a polynomial 
since if 6^ = = bk~i, then so must bk~2 and all the previous coefficients. 
Choose ko so large that for A; > /cq > 4 



ek 1 
l + e(k + 2) - 2 



2 + e{k + 2) 



ek 



< 2 



We choose a specific p large enough so that 



P 



(64) 



(65) 



Since X^fc^fco ^k'^^ is an entire function, it follows that for any /cq, there 
exists a constant M > so that 



P \bk\ < M for k>ko 



(66) 



We redefine M in the relation (j66p to be the least upper-bound for p^\bk\. 
Note that since P{uj) is not a polynomial, M cannot be zero. We now 
introduce 

p'' h 



dk 



M 



and rewrite recursion relation for > /cq + 2 as 

l + ek k-2 
"fc = , , , , p dk-i 



k{2 + efc) ' 2 + ek 
The bound ([M]) translates into 

|c?fc| ^ 1 for k > kQ 



p'd. 



k-2- 



(67) 



(68) 



(69) 



Solving (168]) for d^-i, and shifting the indices fc — > /c + 2, we obtain for 

k > kQ 



dk 



+1 



p(A; + 2) 



2iy l + e{k + 2) 



dk + 



2 + e(fc + 2) 
ekp^ 



dk 



+2 



(70) 



Since M was the least upper bound for p^bk, it follows that there exists 
some k^ > k^ so that 

\dkA>\ (71) 



On using ([MD, ([65]) and 1 4, +2 1 < 1, it follows from ^ that 



> (/c* + 2) 



1 

1-4 



> 1 



(72) 
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which is inconsistent with (I69p . 

We note that this argument not only excludes imaginary eigenvalues, 
but it also shows that except for the trivial translation mode Pq\u}) = 1 
corresponding to A = 0, there are no eigenfunctions regular at all three 
points to = 0, ±1. This follows from simply replacing ly by —iX (or |A|, 
respectively) in the above analysis without necessarily restricting A to be 
imaginary. 

5 Calculation of negative eigenvalues for e > 

We now concentrate on the infinite discrete set of negative eigenvalues An(e) 
which continue the eigenvalues An(l) = —n found in [371 [HH]; here the general 
case of e > is considered while the limit e J, is subject of section [6l 

Observing the parametric dependence of the operator L in (j29p on e and 
A, any solution /5 of the homogeneous equation L{e, \)[3 = is analytic in 
A and e, with the possible exception of e = 0. Therefore the Wronskian of 
any two solutions is an analytic function of (A,e), except at e = 0. Since 
each eigenvalue A,i(e) is determined as a zero of a particular Wronskian, 
it follows that it will change continuously with e, except at e = 0. Since 
eigenvalues {An(e)}„ are all real at e = 1, as e is decreased continuously 
from 1 towards 0, the only way eigenvalues can become complex is through 
collision of erstwhile real eigenvalues, i.e., through the existence of a higher 
order zero of the Wronskian for some e. Such collisions are not observed in 
our numerical calculation, consistent with the fact that all eigenspaces are 
one-dimensional, as is obvious from the behavior of I3^\u}) for a; — > ()30p . 
This suggests that eigenvalues with negative real parts and nonvanishing 
imaginary part are not possible, and therefore they will not be considered 
in the ensuing. 

For the present problem the calculation of A„(e) as zeros of the Wron- 
skian is feasible only for e not too small. Relying on the numerical solution 
of the ODE (|28p . with decreasing e this method rapidly breaks down since 
for |A| <C 1/e, the second independent solution near uj = 1 shows only a 
very weak singularity, cf. Eq. (j3ip . It therefore needs extreme numerical 
precision to determine reliably the solution that is regular at a; = 1. 

To circumvent this problem, we note that in the parameter space spanned 
by (e,A), there exist special points where both independent solutions of 
Eq. (j3ip are regular at u; = 1. These points are found on curves 

A=--m, -<mGN, (73) 

where the general solution near lo = \ can be written as [U] 

/3j)(a;) =ci{uj- iTBiiuj) + C2(^(e) (1 - oj)"^ Bi{uj) \n{uj - 1) + B2{uj)^ . 

(74) 
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Here C1.2 are arbitrary constants, i3i^2(l) = 1 and both -Bi(w) and -62(0;) are 
regular at w = 1. At the zeros eg of A{e), the singularity vanishes. For the 
corresponding A = 1/eQ—m, both independent solutions are regular at u; = 1. 
Therefore solutions regular both at u; = and at u; = 1 can be constructed, 
and A is an eigenvalue for the particular value eo- In our case, these special 
points in the (e, A)-plane can be determined as roots of polynomials in e with 
integer coefficients; and therefore they can be determined with unlimited 
numerical precision. 

We base the calculation on the formal Taylor expansion about uj = 1: 



(75) 



fc=0 



The ODE (j28p yields the recursion relation 



dk 



ek{k - 2) 4„2 - [2A(1 + ek) + 3efc(fc - 1)] 4_i 
2A;(l-e(A + A:)) 



for A: > 1, (76) 



where we take do = 1, d-i = 0. For A as in (j73|) and for k = m, the 
denominator vanishes and for general e the ansatz ()75p breaks down, showing 
that picks up the singular part with A{e) / in ([Till . However, if 

the numerator in Eq. (j76|) vanishes, the solution stays regular. 

In evaluating this condition, it is preferable to rewrite the recursion 
relation (|76|) in terms of 



Pk{e) = e''kl\ll{l-e{X + j))] 4 



(77) 



with A = 1/e — m inserted. This yields 



Pk{e) = e\ik - l)ik - 2){m - k + 1) Pk^2ie) 



k 



1 + e{k - m) + e^- {3k - 2m - 3) 



Pk-i{e) (78) 



with -P-i(e) = and -Po(e) = 1- -Pfc(e) is a polynomial in e of degree 2k, 
with integer coefficients. Any real zero eo of Pm{f) yields an eigenvalue 
A = 1/eo — m. 

Using this approach we determined the 1 1 largest negative eigenvalues for 
values of m ranging from m = 2 to m = 2000, the latter value corresponding 
to e ~ 5 - 10~^. We note that the higher coefficients of Pm{^) with increasing 
m become extremely large, but as function of e, -Pm(e) oscillates around 
zero with an amplitude that for large m becomes extremely small in the 
range of interest < m — 1/e = 0{1). For m = 2000 this amplitude is 
of the order 10"^^'^'^. Nevertheless the zeros of Pmie) can be determined 
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with arbitrary precision since the coefficients of Pm{^) are integers, exactly 
determined from the recursion relation (|78p . Our results are shown in Fig. [1] 
in a double logarithmic plot. For n fixed, |An(e)| is seen to decrease with 
decreasing e. For e | 0, the results suggest the behavior 

A„(e) = -a„ei/2 J^l + 0(e)), (79) 

reflecting the fact that the terms in £e (HHI) that involve second derivatives, 
for e I play the role of singular perturbations. The coefficients On extracted 
from our data are collected in the middle column of Table [TJ In the next 
section we consider the limit e | more closely. 




Figure 1: Double logarithmic plot of the eigenvalues A„(e), n = 1, . . . , 11 
(from bottom to top) as a function of e in the range 5 • 10"^ < e < 1. 
The black dots indicate discrete eigenvalues of the form (j73|) . gray (online: 
red) dots give eigenvalues calculated from the Wronskian. The lines are 
interpolations of the form (|79p . The coefficients a„ determined from these 
lines can be found in Table [H 
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n 


(Xn from Section [5] 


OLyi from Section 16.11 


1 


1 r\ 1 o 1 


1.91 


r) 


4.t30lO 


4.35 


6 


7.oiU7 


7.31 


A 

4 


ID. I loo 


ID. ( Z 


r 

5 


14.5z49 


14.0 


ft 
u 


1 8 f!Q1 fi 

lO.Uc/lU 




7 


23.1907 




8 


27.9992 




9 


33.0985 




10 


38.4732 




11 


44.1103 





Table 1: The coefficients for n = 1, . . . , 11 in the eigenvalue presentation 
Xn{e) = —an e^^"^ (1 + 0{e)) from Eq. ([79]) . The middle column is extrapo- 
lated from Section [5] and Figure [H The right column is from section 16.11 

6 The eigenvalues for e | 
6.1 Eigenvalues with low indices 

According to Eq. (|31|) . the general solution of the ODE (|28|) at = 1 
develops a singularity that becomes arbitrarily weak if e tends to zero. As 
a consequence, representing the special solution (3reg{^) regular at w = 
as a power series in e and imposing regularity at a; = 1 does not put any 
constraint on the eigenvalue A. To show this we write firegi^) in the form 

oo 

Preg{co)=e^^J2''M^)^ (80) 
3=0 

imposing the normalization condition 

PregiO) = 1, 

i.e. 

/o(0) = l, /,(0) = 0, i>l. (81) 
Substituting this ansatz into Eq. ([28]) we find 

Ulo) = 1 (82) 

fj(u!) = -1(1 - UJ^) UJ d^fj-i{uj) + XuJ^fj-l{uj) 

+^ j dx {l + \x- 3x^ + \x^) fj-iix) for j>l. (83) 
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Evidently fjiuj) is a polynomial in uj of degree 4j. Thus (3reg{^), Eq. (fSOl) . 
evaluated to arbitrary finite order in e, is analytic at tj = 1 for arbitrary 
A. This breakdown of an approach based on a simple expansion in the 
regularizing parameter e is a well known feature of such moving boundary 
problems, both in steady state (see for instance [8]) and linear stability 
analysis (see for instance |12]). 

For later use we note the expression for the low order terms: 



f^{u:) = -^(l + 2a;2_3^4)^Ac^(_i^3^2^5^4) 
4 o 



X^oj^ 23 , 15 A AV 



4 4/ 4 



We also note the singular solution of Eq. ()28p . resulting from a WKB-analysis 

/3.m,(c^) = (l - ^) (tT^) ' + • ^^^^ 

This solution picks up the singularities both at w = and at uj = ±1. 

The above analysis demonstrates that for e | the eigenvalues A can be 
determined only by an analysis going beyond all orders in e. Indeed, taking 
into account the behavior A ~ eV2 sug gested by section El it is evident that 
the individual terms of the asymptotic expansion (jSOp become of order 1 for 
\uj\ ~ e~^/^ and consequently the expansion becomes invalid. Since, however, 
the eigenfunctions fi^^"^ {uj) by definition are holomorphic in the right half uj 
plane, this suggests to analyze the region ~ e""*^/^ more closely. 

We therefore introduce the scaled variable 

7 = (86) 

and the notation 

A = (87) 

/3j)(a;(7)) = oil). (88) 

With these substitutions, the operator L(e, A) from ()29p takes the form 

L(e,A) = — 7 (Lo(a)-eLi(a)j, (89) 



Lo(«) = 9^+(27-^)5,-^, 
Li(a) = 7292 + (^_2a)a^ + ^. (91) 



(90) 

2a 

7 
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To leading order in e we have to discuss the ODE 



^o(«)5(7)=0. (92) 
An asymptotic expansion in powers of ^ yields the solution 

s.(T)~i-- + j^-j^(i + Yj+o(:^). m 

The second solution, found by balancing the terms 9^ and (27— ^) 3-^, takes 
the form 

52(7) {\^o(^ . (94) 

It is easily seen that in the range 

e~^/^ » I7I » 1, i.e., 1 < < e~^/^ , 

^'1(7) matches the regular solution /3reg{^^) given in ([50]) and ([5^ . whereas 
^2(7) matches the singular solution figingi'^) given in (fSSl) . 

Evidently 52(7) dominates over 51(7) in the two sectors ^ < | arg 7I < |7r 
for both signs of arg 7. We search for eigenfunctions /3(a;) matching (71(7) in 
the complete range | arg7| < ^. Now it is well known that starting deep in 
the region ^ < arg 7 < ^ with initial conditions taken from the asymptotic 
expansion (|93p and integrating Eq. (|92p down into the region — § < arg 7 < 
— ^, we in general will pick up a dominant contribution proportional to 
5(2(7). If, however, 51(7) G M for 7 S M, the Schwarz reflection principle 
guarantees the absence of such a contribution. Equivalently, we may state 
that by definition the eigenfunctions f3{ijj) are real for uj € (0, 1), and that 
the absence of the singularity induced cut for a; > 1 implies /3(<-o') G M for 

a; > 0. Thus the eigenvalues A^^ = —an(^^'^{^ + 0{e)) are selected by 
imposing the constraint (7(7) G M for 7 > 0. 

To evaluate this criterion we started at points 70 = b+iL, b = 0.1, 0.2, 0.3, 
L = 20, 25, 30 and integrated Eq. (j92p down to the real axis along lines 
b = const. The results for the first 5 positive zeros a„ of Im[g(6)] are given 
in the right column of Table [TJ They clearly are consistent with the asymp- 
totic results found in subsection [5] and presented in the middle column of 
the same table. We did check that the quoted values to the precision given 
are insensitive to the starting value 70 as long as L is chosen sufficiently 
large and 6 > is not too large or too small. The former causes Im[(7(6)] 
to be rather small and zeros are harder to detect, while the latter causes 
numerical inaccuracies due to the singularity of LQ{a) at 7 = 0. Fig. [2] 
shows that with different b, while Im[(7(6)] itself is different, the zeros match 
as expected from theory. 
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6.2 Eigenvalues with large indices n 



Both the approach of the last section and that of Sect. [5] within reasonable 
numerical effort yield the asymptotic coefficients an only for the first few 
eigenvalues, On = 0{1). However, in the complementary region ^ 1, a 
fully analytical analysis is possible, similar to that performed earlier in the 
surface tension selection problem for steady Hele-Shaw fingers [HI \TT\ I13j . 
We first introduce a Liouville transformation to eliminate the first order 
derivative in the operator Lo(a). With 



ail) 



7I/2 e-^V2 



Eq. ([921) takes the form 



2a 



+ 



472 



h = 0. 



We now rescale 7 according to 

7 = a^/^x, Klix)) = q{x) 



to find 



a 



4/3 



q = 0. 



For a — > 00, a WKB- analysis yields the asymptotic relation 
q{x) ~ Ci qi{x) + C2 92 (x), 



91,2 (x) 



Q-'/\x) exp 



^0 



(95) 

(96) 

(97) 
(98) 

(99) 
(100) 
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where 

Q(x)=x' + -. (101) 

X 

92 (x) dominates in the sectors ^ < | argxl < f , but is subdominant on the 
real axis. Recalhng Eqs. (f95|) and ([971) . it is easily seen that for large 
qi{x) yields 51(7), Eq. ([93]) . whereas (?2(x) yields 52(7)) Eq- (jMj) . Since the 
eigenfunctions in the limit e J, in all the right half plane Re [7] > must 
reduce to g\ (7) , this implies that Ci is real whereas C2 has to vanish in both 
sectors f < | argxl < f • 

Now the WKB-analysis breaks down at the turning points Xs defined by 
Q{Xs) = 0, which in the relevant region Re[x] > yields the two solutions 

= 2^/3 e^'^/'\ (102) 

As is well known in asymptotics [58] the coefficients in the asymptotic rela- 
tion (I99p may jump when x crosses a Stokes line emerging from a turning 
point. Thus a vanishing C2 in the sector | < argx < f does not im- 
ply that C2 vanishes on the positive real well. Since, however, the 
eigenfunctions are real on the real axis, we get the relation 

Im[C2] =0 for X G IR+ (103) 

as a necessary condition fixing the eigenvalues. To evaluate this condition 
we determine the jump of C2 by analyzing the neighborhood of the turning 
point Xs = 2i/3e+W3. 

Expanding Q(x) about Xs, 

Q(x) = Q'iXs) {x-Xs) + 0{x-Xs)\ (104) 

Q'iXs) = (105) 

and rescaling x~ Xs according to 



e = e-^'^^'^ix-Xs), (106) 

q{xm = HO, (107) 

we find that Eq. (j98p to leading order in a reduces to the Airy equation 

(5|-e)m(e) = 0. (108) 

The phase in Eq. ()106p has been chosen such that for ^ > we enter the 
region where C2 = in the asymptotic relation (j99p . This analysis is valid 
for Ix ~ Xs| ^ a~'^/^ , and using Eq. (I104p . we find for small |x — Xs\ along 
the line ^ > 



/(x) ~ a r'/' 



exp 



2^3/2 

3^ 



(109) 
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where the constant a is given by 

1/6 



a = Ci 



22/3 

~6~ 



e-^2V9 exp 



a 



2/3 



Jo 



(110) 



This result is valid for a^/^ ^ ^ 1, where it must match the solution of 
the Airy equation (jlOSp . This yields 



m 



(0 = 2y^aAi(0. 



(Ill) 



Continuing this result clockwise around the turning point to negative ^, we 
find 



m(e) ~ a ler^/^ e^"/^ (exp 



I exp - i|4r'' ) for — ^ ^ 1. 

(112) 

For a^/^ > > 1, this again has to match with Eq. (|99|) , which in this 
region reduces to 



exp 



3/2 



(113) 



C2 
+ ^ exp 



2a2/3 / dx'Vol 



X' 



exp 



3/2 



We thus find 



C2 = —iCi exp 



2a2/3 rdx'VOM 
Jo 



(114) 



Since both Ci and C2 have to be real, this leads to the eigenvalue condition 



Im 



2a2/3 / dx'VOM 







n + ^ I vr, n G 



(115) 



The integral is easily evaluated to yield 



dx'VOix^ 







2^/3 e-/6 / ds^ --s^ 



22/3 



V^r(i; 



8 r(|) 



(\/3 + i). 



(116) 



Eq. (|115p yields the final result 



ar, 



3/2 



3/2 

1.13254... l^^+2) forn G Nq. 



(117) 
(118) 
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Figure 3: from Eq. ()119p as a function of n, together with the inter- 
polating formula p2U|) . The plot demonstrates the accuracy of the asymp- 
tote (I117p even for small n. 



By construction this result is valid in the limit oo, i.e., for n oo. 

However, it turns out to be a good approximation for small n as well. This 
is illustrated in Fig. [3] where we plot 



with an, n = l,...,ll taken from the second line of Tabled) Even for n = 2 
the asymptotic result (jllTp differs from the true eigenvalue by less than 3 
%. The interpolating curve included in Fig. [3] is given by 



indicating that the asymptotic limit (jllZp is rapidly approached. 

7 Behavior of low order eigenf unctions for |<x;| < 1 

As is evident from Eq. (|3ip . for w — > — 1 the eigenfunctions are singular, 



For small e and |A| <^ 1/e, this singularity will show up only in derivatives 
of high order, and the Taylor expansion (I49p will converge in the whole 
physical domain < 1. We therefore can use this expansion together 




■3/2 



(119) 




(120) 
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Figure 4: The ratio T = (3^^^ (w) je^^ 
for e = 0.1, 0.05, 0.025 as indicated, a) Absolute value \T\ and b) argument 
Arg(.F)/7r as a function of a/vr. 



'^^^ evaluated on the unit circle a; = e*" 



with the recursion relation (j50p . (j5ip to evaluate pf_ \uj) for the first few 
eigenvalues. Outside some neighborhood of a; = —1 these functions are 
expected to govern the large time behavior of analytic perturbations. 

^From the discussion of P^egi^) in Subsect. 16.11 we expect that in the 
physical domain |tij| < 1 the eigenfunctions for small e are well approxi- 
mated by simple exponentials. This indeed is true as illustrated in Fig. HI 
Panel mia) shows the absolute value \p{^^{uj) e~'^^'^\ evaluated on the unit 
circle lo = for a decreasing set of values of e. Evidently this ratio tends 
to 1 and is not very different from 1 even for e = 0.1. Panel H^b) shows 
the corresponding phase difference which is found to be very small and to 
tend to zero with decreasing e. All curves in Fig. H] are well represented by 
the expansion ([HOjl evaluated up to order e^. As expected, for fixed e the 
approximation 

/?J)(e*") ~ exp(Acosa + iA sin a) (121) 

becomes worse with increasing A. This is illustrated in Fig. [5] for e = 0.1, 
Ai = -0.5462, Ae = -4.6086. Panel a) compares Re ^In /3|'^(e*°)) /A 
with cos a. In the neighborhood of a = vr the deviation from cos a increases 
with increasing A, which is a consequence of the singularity of w = —1. But 
outside that range the approximation ()12ip is quite accurate even for ■ 

Panel mb) shows the corresponding results for Im ^In Px^^{e^°')^ /A. 

As pointed out in section [231 we expect that the initial condition P{uj, 0) 
in some neighborhood of u; = 1 can be expanded in terms of eigenfunctions 
as 

oo 

/3(..,0) = 5^5n/32M, (122) 

fc=0 

and that the domain of convergence of this expansion increases with time 
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Figure 5: Real part (a) and imaginary part (b) of the function Q 



In 



/3|'^^(e*°) /A for e = 0.1 for Ai (black solid line) and Ag (blue or grey 

solid line) as a function of a/vr. The dashed lines show cos a in (a) or sin a 
in (b) for comparison. 



to ultimately cover U^^ \ {—1}. For small e and not too large |A|, (3^ {uj) 
is well approximated by e^^ , and all eigenvalues vanish for e — > 0. Thus 
for fixed uj and e — > 0, all eigenfunctions tend to 1. This suggests that for 
small e many terms must contribute non-negligibly to the expansion of a 
generic initial condition. Furthermore the phase of the coefficients Qn must 
vary strongly. With decreasing e it will need more and larger coefficients 
to describe a non-trivial general initial condition. These large coefficients 
imply large transient growth, because in the representation 

oo 

/3(.;,t) = J]5n/3g(^)e""*, (123) 

fc=0 

the factor e'*'"* is substantially different for each n when t = 0(e~^/^). The 
large transient growth restricts the validity of linear stability analysis to a 
small ball of initial conditions. This is conffi'med in part II in the study of 
the initial value problem. 

8 Conclusion and outlook 

In the introduction we raised the question: Can a kinetic undercooling 
boundary condition regularize the evolution of a curved interface in a Lapla- 
cian growth model? This problem is motivated by the physics of streamer 
discharges which determine the early evolution of sparks and lightning. We 
here gave a first answer by analyzing the spectrum of linear perturbations 
of uniformly translating circles. We proved that for all e > the spectrum 
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is discrete and that all eigenvalues A„,, except for Aq = 0, have a negative 
real part. Thus any infinitesimal perturbation tends to a constant expo- 
nentially in time for large time, and asymptotically the circular shape is 
recovered. For arbitrary regularization parameter e > we found an infi- 
nite set of negative real eigenvalues A„(e); they smoothly continue the exact 
result A„(l) = —n, n G N, found previously for e = 1 [38]. 

In formulating the eigenvalue problem we had to allow for a singularity 
of the eigenfunctions at the point a; = — 1 at the back of the circle. Therefore 
an expansion of a regular initial perturbation in terms of eigenfunctions must 
break down in a neighborhood of = — 1. For e = 1, it was found [38j that 
the size of this neighborhood with increasing time r decreases as e~'^, and 
that all the structure of the initial perturbation is convected into this region. 
In part II [l5] of this series of papers we will argue that such behavior is 
found for all e > 0. 

Our results suggest that in the framework of linear perturbation theory 
the evolution of a curved front can be regularized by a kinetic undercool- 
ing boundary condition. Linear perturbation theory must break down in 
the limit e | 0. In our results this shows up in the asymptotic behavior 
of the eigenvalues A„(e) = —an e^''^, where a„ ~ const n^^'^ for large n, 
and in the associated behavior I3^\uj) — > e'^"'^ — > 1 of the eigenfunctions. 
This indicates that for small e, the eigenfunction expansion of even a very 
smooth initial perturbation of small but finite amplitude will contain many 
terms with large coefficients. Initially these terms almost compensate, but 
the balance is destroyed by the temporal evolution. Generically this will 
lead to a strong transient growth of the perturbation which may drive the 
evolution into the regime where nonlinear effects have to be included. Ex- 
amples supporting this scenario will be given in part II }45] . For e ^ 1 this 
mechanism can lead to an instability of the circular shape also against quite 
small and smooth perturbations. 
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A Appendix A: Asymptotics of Taylor series co- 
efficients and relation to singularity at u = Us 



Here we prove the following standard results for the sake of completeness. 
Theorem: Assume f{uj) is analytic in \lu\ < R, except for a singularity at 
a; = / (with \uJs\ < R), where 

f (til) = f situ) {UJ - UJs)°' + faiuj) as UJ^tHs. 

Here a is not a positive integer, and fg and fa are locally analytic at lo = uJs 
with fs{i^s) ^ C / 0. Then the Taylor series coefficient bk in the represen- 
tation 

oo 

f{u;) = Y,hu^'' 

k=0 

has the leading order asymptotic behavior 

h ^sm(7raK fc.+i ■ 



Proof: We will first assume q > — 1. We recall the contour integral repre- 
sentation of 6fc 



=5 ^ 



where 6 is chosen sufficiently small so that \uj\ =5 contains no singular- 
ity of /(w) and \ujs\ + S < R. We deform the contour into f^^^ + /^^ + f^^ 
as shown in Figure [6j The contribution from J^^ is easily bounded by 

M [\ujs\ + f^]"^, where M = sup^^^^^^^^^g \ f{uj)\. Now consider the contribu- 
tion from J^^ + J^^ . It is convenient to introduce a change of variable 



C = log ■ 



Then it is readily checked that we have 

r.<5i 

I UJ ^ I / 

+ 



where 5i = log 
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(124) 



Noting that the contribution from fa cancels out 

/,(u;(C))e"^<dC- 



We note that in the neighborhood of ^ = 0, 
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So, using Watson's Lemma, we obtain 

C 



hi 



■k+a 



[1 



k 



r(i + Q) 



from which the Lemma follows since the contribution from J^^ is evidently 
exponentially small relatively for large k. 

If a < —1, we consider n-th iterated integral Inf, where [/i/](u;) = 
f{uj')dio', = Ii (hf) and so on. Then, it is clear that the singularity 
of /„/ at to = LVg will be of the type {to — uJs)"'^"', and we can arrange 
a + n > —1. The argument above can then be repeated for power series 
coefficients of /„/. The Theorem follows on differentiating n-times the power 
series of /„/. 

Remark The previous theorem remains valid even when a is a negative 
integer, provided the product of sin(7ra)r(l + a) is replaced by its finite 
nonzero limit lima_>_.„. The validity of the result is easily checked from 
Taylor expansion of ^^^^ ^ and its derivatives in terms of a geometric series 
and its derivatives. Further, since the asymptotic result relies on Watson's 
Lemma, the condition on analyticity of fs at u = ujg can be weakened to 
fs{u;) = C + 0{{co-LOsf) forP>0. 




Cq:={co: |co|=|cOg| + 6 , v < arg co< v+ln } 

Figure 6: Contour ^^1=5 deformed to /^^^ + /^^ + J^^ for evaluation of bk- 



Theorem: If the k-th Taylor series coefficient b^ of an analytic function f 
at u! = satisfies the following 

b, ~ -^e-" sin(™)^7'=+-%±i^ (1 + 0(1/^^)) (125) 
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for non-integral a, then 



f{u;)r.C{LO-iOsr + fa{co), (126) 

where fa{uj) is regular atuj = uJs- The same conclusion is valid for a = —n, a 
negative integer, provided r(l+a) sin(7ra) is replaced by its limit as a ^ —n, 
On the otherhand, if 

bk^ -Coj-''k-'{l + 0{l/k)), 

then 

/(a;) ~ Clog (l-;^)- 

Proof. We will assume for now that a E (—1,0). Then using previous 
Theorem to determine Taylor series coefficient of C(w — Wg)", it follows on 
subtraction that 

oo 

g{uj) ^ f{uj) - C{uj - LO.r = Yl ^k^" 

will have Taylor series coefficient bk ~ Const. u>J^k~°'~'^ for large k implying 
that the series is absolutely convergent on = \iOs\. In particular, g is 
continuous at u; = uJs, which proves the theorem for a € (—1,0). If a takes 
on other ranges of values, we obtain this result by either n-times iterative 
integration from the origin of the power series for /, or differenting it n times 
so as to ensure that for non-integral a, a + n £ ( — 1, 0) or a — n G ( — 1, 0). 
The result quoted in the Theorem follows by noting that the asymptotics 
is differentiable as it is valid for a; — in a complex sector. The second 
result follows from noting the explicit Taylor expansion of log ^1 — and 
noting that the difference has a series that is absolutely summable and hence 
the remainder is continuous at lu = ujg- By using explicit derivatives or 
integrals of Geometric series, the conclusion (I126p holds for for negative 
integer q = — n as well. 
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